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INTRODUCTION 


V 


Theoretical  models  and  computer  simulations  of  the  interaction 
between  narrow-band,  whistler-mode  waves  and  energetic  electrons  using 
parameters  suitable  to  the  magnetosphere  show  that  the  resonant  electrons 
are  phase-bunched  by  the  waves  [e.g.  Brinca,  1972;  Helliwell  and  Crystal, 

1973;  Nunn,  1974]  . In  a monochromatic  wave,  the  motion  of  an  electron 
is  determined  by  the  force  acting  on  the  particle  from  the  electric  and 
magnetic  fields  of  the  wave,  A feedback  loop  is  established  between  the 
radiation  emitted  by  the  phase-bunched  electrons  and  the  resonant  electrons 
entering  the  interaction  region.  The  dominant  force  is  the  magnetic  force 
which  causes  a longitudinal  drift  of  the  electrons,  i.e.  a drift  along  the 
geomagnetic  field  line. 

If  9 is  the  angle  between  the  electron  velocity  component  v^  perpendicular 
to  the  magnetic  field  and  the  magr  field  component  Bw  of  the  wave,  the 
longitudinal  drift  produces  an  effective  phase  rotation  between  vA  and  B^. 

This  phase  rotation  changes  the  force  responsible  for  the  longitudinal  drift. 

The  "equation  of  motion"  for  the  angle  0 is  [ Helliwell  and  Crystal,  1973] 


d20/dt2  + (-3—)  k v.  B sin  6 = 0 
\ m / i w 


! I 


(i) 


This  is  the  same  form  as  the  equation  for  the  anharmonic  motion  of  a ffmple 
pendulum. 

If  the  phase  of  the  transmitted  wave  is  changed  abruptly,  the  solution 
to  (1)  requires  a new  initial  value  for  9 established  by  the  phase  change  of 
the  wave. 
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Controlled  experiments  to  observe  the  effect  of  phase  changes  of 
a whistler-mode  wave  in  wave-particle  interactions  in  the  magnetosphere 
were  undertaken  in  conjunction  with  very -low -frequency  (VLF)  wave  trans- 
mission experiments  conducted  at  Fort  Richardson,  Alaska  in  the  summer 
of  1973,  These  experiments  were  performed  with  the  transportable  very- 
low-frequency  (TVL.F)  transmitter  facility  [Koong  and  Dazey,  1974], 
Between  September  3 and  September  14  phase  reversal  programs 
were  transmitted  daily  at  frequencies  of  6.6  kHz  and  13.275  kHz.  At 
13.275  kHz  the  TVLF  transmitter  radiates  approximately  1000  W.  The 
primary  ground-based  receiver  in  the  conjugate  area  was  located  at 
Dunedin,  New  Zealand  (L  - 3).  At  that  site,  the  amplitude  of  the  earth- 
ionosphere  waveguide  signal  at  13.275  kHz  was  typically  10  to  20  pV/m. 

For  cw  type  transmissions  this  made  it  difficult  to  separate  the  magneto- 
spheric  signal  from  the  waveguide  signal  which  was  typically  a few  pV/m. 

At  6.6  kHz  the  radiated  power  was  100  W.  With  the  increased  attenuation  at 
the  lower  frequency,  the  waveguide  signal  was  generally  below  the  receiver 
sensitivity.  The  data  discussed  below  were  obtained  during  transmissions 
at  6.6  kHz. 
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The  TVL.F  transmitter  is  a linear.  Class  AB,  power  amplifier 

I * 

driven  by  a low-voltage,  highly -stable,  frequency  synthesizer.  A 


block  diagram  of  the  transmitter  and  control  circuit  including  the  phase 
reversal  switch  is  shown  in  Figure  1.  Frequency  stability  is  provided 
by  a 5 MHz  quartz  crystal  oscillator  in  the  frequency  synthesizer  which 
generates  the  6.6  kHz  signal  for  the  transmitter.  The  stability  of  the 

g 

oscillator  is  ± 1 part  in  10  per  day.  A 1 MHz  signal  derived  from  the 
5 MHz  oscillator  also  drives  a time -code  generator.  The  time-code 
generator  is  set  within  ±2  msec  of  Universal  Time  by  synchronization 
of  a 1 pulse-per -second  output  with  the  second  marks  transmitted  by  the 
Time  Standard  Station  WWVH.  Proper  account  was  taken  of  the  propagation 
delay  from  WWVH  to  the  TVLF  transmitter  site. 

An  identical  frequency  synthesizer  at  the  receiver  site  in  New 
Zealand  was  synchronized  with  the  Australian  Time  Standard  Station  VNG. 
The  time  code  recorded  at  the  receiver  site  was  also  maintained  within 
± 2 msec  of  Universal  Time  and  the  overall  system  timing  accuracy  was 
± 4 msec. 

At  periodic  intervals  phase  inversion  (180°  phase  shifts)  of  the 
driving  voltage  is  achieved  by  a COSMOS  integrated  circuit  reversing 
switch  that  is  keyed  by  a programmer  driven  by  the  time -code  generator. 
The  switch  accomplishes  the  phase  reversal  in  about  1 psec.  It  can  be 
keyed  32,  16,  8,  4,  2,  1,  0.5,  0.2,  0.  1,  0.0333...,  or  0.01666...  times 
per  second.  The  signal  out  of  the  switch  passes  through  a unity  gain 
differential  amplifier,  an  audio  amplifier  and  finally  provides  the  input 
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signal  to  the  Intermediate  Voltage  Amplifier  of  the  transmitter. 

The  main  instrumental  effect  attributed  to  the  transmitter  is  the 
rise  time  of  the  tuning  circuit.  The  antenna  system  is  essentially  a 
series  RL.C  circuit  fed  by  a matching  transformer.  The  circuit  is 
shown  in  Figure  1.  The  current  to  the  antenna  is  measured  by  a calibrated 
current  probe  on  the  output  side  of  the  matching  transformer.  The  Q of 
the  circuit  was  periodically  measured  by  sweeping  the  transmitter  fre- 
quency about  the  center  frequency  and  measuring  the  bandwidth  at  0.707 
of  the  maximum  current.  About  25  minutes  after  the  phase  reversal 
transmissions  discussed  below,  the  Q was  measured  to  be  68  which 
corresponds  with  a current  rise  time  of  3.  3 msec. 
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OBSERVATIONS 


On  September  6,  1973  strong  whistler-mode  signals  from  the 
transmitter  were  detected  at  6.6  kHz.  During  a 20-sec  period  centered 
at  21:43:30  Universal  Time,  strong  signals  were  detected  while  a phase 
reversal  program  was  being  transmitted.  The  signal  amplitude  averaged 
over  a 10-sec  period  was  18  pV/m.  The  noise  in  a 50  Hz  bandwidth  at 
6500  Hz  was  3 pV/m.  During  this  time  period,  phase  reversals  were  keyed 
at  0.5  sec  intervals.  Figure  2 shows  a record  of  the  signal  amplitude  as  a 
function  of  time.  During  that  time  period  the  amplitude  of  the  received 
signal  gradually  increased  10  db.then  returned  to  its  previous  level. 

In  Figure  2 it  can  be  seen  that  the  signal  amplitude  periodically 
decreases  in  amplitude  and  then  recovers.  This  occurs  every  0.5  sec 
which  is  the  time  interval  at  which  the  180°  phase  shifts  were  being  applied 
to  the  transmitted  signal. 

In  Figure  3 a short  section  of  the  amplitude  record  is  reproduced  with 
higher  time  resolution  together  with  a measurement  of  the  phase  of  the 
received  signal.  The  phase  measurement  shows  an  abrupt  phase  dis- 
continuity of  180°  each  time  the  amplitude  goes  through  a minimum.  These 
measurements  demonstrate  that  phase  coherence  is  maintained  throughout 
the  magnetospheric  propagation  path. 

In  Figure  3 a sinusoidal  amplitude  modulation  is  also  apparent 
particularly  in  the  signal  at  12:43:34  UT.  Most  of  the  record  is  contaminated 
by  impulsive  sferics.  However,  no  sferics  were  apparent  in  the  data  during 
the  0.5  sec  time  period  between  the  phase  reversals  which  straddle  12:43:34  UT 
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Amplitude  of  the  6.6  kHz  signal  received  at  Dunedin,  New 
Zealand  as  a function  of  time  during  a period  when  the 
transmitter  was  sending  a program  in  which  the  phase  of  the 
voltage  driving  the  power  amplifier  was  reversed  at  0.5  sec 
intervals.  The  amplitude  scale  is  linear  in  magnetic  field 
intensity. 
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The  decreases  and  increases  in  the  amplitude  of  the  received 
signal  at  each  phase  reversal  occur  approximately  exponentially  with 
equal  time  constants,  t , An  example  of  the  signal  amplitude  plotted  as 
a function  of  time  during  a period  of  decreasing  amplitude  is  shown  in 
Figure  4,  The  time  constant  obtained  from  the  slope  of  the  line  is 
33  msec. 

In  order  to  obtain  a reasonable  signal-to-noise  ratio,  a narrowband 
filter  was  required  in  the  analysis  instrumentation.  Figure  5 shows  the 
result  of  passing  a short  segment  of  the  data  through  three  filters  with 
bandwidths  of  80  Hz,  240  Hz,  and  800  Hz.  The  decrease  in  signal  amplitude 
at  0.5  sec  intervals  is  clearly  seen  (at  0.13  sec,  0.63  sec  and  1.13  sec)  with  all 
three  filters.  The  total  time  for  the  signal  to  decay  and  recover  is  about  0.  1 sec. 

It  is  interesting  to  note  that  neither  self-sustaining  oscillations  at 
the  previous  phase,  nor  triggered  emissions  were  observed  at  the  time  of 
a phase  reversal. 

Two  distinctive  features  of  the  signal  which  will  be  discussed  are 
the  rise  and  fall  times  of  the  signal  amplitude  at  each  phase  reversal  and 
the  amplitude  modulation  between  phase  reversals.  The  time  constants 
and  periods  observed  are  characteristic  of  the  interaction  times  of  electrons 
with  whistler-mode  waves.  They  are  significantly  longer  than  those 
associated  with  the  bandwidths  of  the  antenna  tuning  circuit  or  the  signal 
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Amplitude  of  the  received  signal  as  a function  of  time  during 
a signal  decrease  following  a phase  reversal. 





DISCUSSION 

Interaction  Bandwidth.  The  decrease  and  subsequent  increase 
of  the  signal  amplitude  at  each  phase  reversal  is  analogous  to  the  response 
of  a resonant  RLC  electric  circuit  when  the  phase  of  an  applied  sinusoidal 
voltage  at  the  resonant  frequency,  w*  = 1 / J LC,  is  reversed.  In  the 

]jt 

electrical  circuit  the  decay  time  r for  the  amplitude  is  related  to  the  Q 

$ 

of  the  circuit  and  to  the  bandwidth,  A « /2ir  by 

Q = w*t*/2  = w*  / A w*  (2) 

This  bandwidth  arises  from  the  dissipation  of  the  energy  by  ohmic 

heating  in  the  resistive  components  of  the  circuit. 

In  the  VLF  phase  reversal  experiment,  the  decay  time  for  the  wave 

amplitude  following  a phase  reversal  is  measured  to  be  33  msec.  For 

this  time  constant,  the  effective  Q obtained  from  (2)  is  684  and  the  bandwidth, 
$ 

Aw  /2ir,  is  9.6  Hz.  Several  "dissipation"  processes  influence  the  effective 
bandwidth  in  the  resonant  wave-particle  interaction.  Among  them  are 
Landau  damping,  phase  mixing  in  the  non-uniform  magnetic  field,  and 
radiation. 

Schulz  [1972,  19743  has  derived  the  following  expression  for  the 
intrinsic  bandwidth  of  cyclotron  resonance  in  a non-uniform  magnetic 
field 


Aw/2ir  = | (w/16ir)  | 1 /3[l  - v„/vg)]"2/3 


(3) 
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where 


«• 

(*) 


MW  | 

“37“  > 


(4) 


where  s is  the  arc  length  along  the  field  line  and  w is  evaluated  at  the 

_•*  A 

equator  (s  = o),  vg  = v • B,  and  vg  * dw/dk, . 

For  the  plasma  parameters  appropriate  to  the  phase  reversal 
experiment  (see  below)  the  intrinsic  bandwidth  is  calculated  to  be  26  Hz. 

One  must  be  cautious  in  comparing  this  with  the  observed  value. 

The  intrinsic  bandwidth  calculated  by  Schulz  [1972,  1974}  includes  only 
the  effect  of  the  non-uniform  magnetic  field  on  the  bandwidth.  Other 
effects  such  as  wave  growth  (which  would  tend  to  increase  Aw)  and  phase 
trapping  (which  would  tend  to  decrease  Aw)  could  alter  the  actual  bandwidth 
from  this  "intrinsic"  bandwidth. 
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Amplitude  Modulation.  Amplitude  modulated  whistler-mode  signals 
from  VLF  transmitters  have  previously  been  reported  by  Likhter  et  al. 

[l97l]  and  Bell  and  Helliwell  [l97lj.  This  modulation  is  believed  to 
occur  at  approximately  the  phase  bunching  frequency  of  the  electrons  which 
are  undergoing  pendulum-type  motion  in  the  potential  well  of  the  wave 
[ Brinca,  1972;  Nunn,  1974,  1975].  However,  Helliwell  and  Crystal  [l975] 
point  out  that  the  amplitude  pulsations  are  a fundamental  result  of  the 
feedback  in  their  model  and  not  a consequence  of  the  oscillatory  electron 
motion  due  to  the  phase  bunching. 

The  modulation  period  scaled  from  the  TVLF  experiment  data  between 
12:43:33.65  UT  and  12:43:34.15  UT  shown  in  Figure  3 is  greater  than  the 
actual  period  because  the  signal  was  passed  through  a 10  Hz  filter  to 
produce  that  display.  The  modulation  period  was  determined  by  measuring 
the  frequency  offset  of  the  sideband  responsible  for  the  modulation.  For 
this  purpose,  maximum  entropy  spectrum  analysis  was  applied  to  a 250 
msec  interval  of  the  data.  The  beginning  of  the  interval  was  12:43:33.69  UT. 
The  data  were  sampled  at  25,  000  samples  per  second.  The  program 
developed  by  Ulrych  and  Bishop  [1975}  to  evaluate  the  Burg  estimates  of  the 
autoregressive  coefficients  was  used  for  the  spectral  determination 

Since  the  maximum  entropy  spectrum  is  an  optimally  smoothed  one, 
Ulrych  and  Bishop  [1975]  conclude  that  the  resolution  of  the  Burg  maximum 
entropy  method  spectrum  is  almost  twice  that  of  a periodogram  spectrum 
[Radoski  et.al. , 1975], 

The  power  spectral  density  for  the  TVLF  data  using  1000  coefficients 
is  plotted  in  Figure  6.  The  points  are  plotted  at  one  Hertz  intervals.  The 
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dominant  peak  occurs  at  the  TVLF  carrier  frequency  with  a maximum  at 
6599  Hr.  The  offset  from  6600  Hz  is  attributed  to  an  imperfect  lockup  with 
the  reference  fA-equency  during  digitization.  An  upper  sideband  with  a maximum 
at  6638  Hz  is  also  present  in  the  spectrum  [ Brinca,  1972]. 

The  difference  frequency  between  the  carrier  and  the  sideband  is 
39  Hz.  The  corresponding  modulation  period  would  be  26  msec.  The 
total  power  in  the  carrier  obtained  by  integrating  the  spectral  density  over 
the  peak  is  6.2  times  the  total  power  in  the  sideband. 


Amplitude  Saturation.  Ossakow  et  al.  [1972]  have  presented  the  results 
of  computer  simulation  studies  of  whistler  instabilities  in  anisotropic, 
collision-less  plasmas.  They  argue  that  wave  growth  saturates  when  the 
bounce  frequency  of  a particle  in  the  potential  minimum  of  the  wave  becomes 
of  the  order  of  the  linear  growth  rate,  w - X,  where  w is  the  particle  bounce 
(or  "trapping")  frequency  given  by  u = 2tt/t  and  y is  the  linear  growth  rate. 
The  computer  simulation  study  of  Ossakow  et  al.  [l972]  yielded  the 
saturation  condition 


where,  for  bi-Maxwellian  distribution  functions,  a is  a number  between  0.25 

and  0.6,  and, for  two  cases  with  Maxwellian  distribution  functions  with  a 

loss  cone,  a is  0.2  and  0.25.  For  a case  representative  of  the  magnetospheric 

plasma  [Ossakow  et  al. , 1973]  a = 0.2.  For  the  computer  simulations,  they 

define  a discreteness  parameter  d = y /k  A V where  k is  a quantized  wave 
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number,  and  AV  is  the  interval  between  two  discrete  adjacent 
res 

resonant  velocities.  When  d =*  1,  then  a particle  is  only  in  resonance  with 
one  wave  at  a time.  This  is  the  condition  required  to  simulate  an  inter- 
action between  the  electrons  and  the  coherent  signal  transmitted  by  the 
TVLF  transmitter.  For  one  of  the  bi-Maxwellian  cases,  d ~ 0.  1 and 
a = 0.3;  for  the  magnetospheric  plasma  case  d ~ 0.2  and  a = 0.2. 

If  we  identify  the  time  constant  of  the  increase  in  signal  amplitude 
following  a phase  reversal  with  the  linear  growth  rate  of  the  wave,  we 
have  a measured  value  of  y = 1/  t ~ 30  sec  * (recalling  that  the  time  constants 
of  the  signal  increases  and  signal  decreases  are  the  same).  This  growth 
rate  is  significantly  higher  than  that  reported  by  Helliwell  and  Katsufrakis 
[19743.  In  comparing  our  results  with  those  shown  in  Fig.  4 of  Helliwell 
and  Katsufrakis  [19743  we  also  note  that  the  time  to  reach  saturation  is 
three  to  four  times  longer  in  their  data  than  in  the  current  experiment.  The 
two  experiments  are  not  entirely  comparable  because  the  particle  phases, 
with  respect  to  the  magnetic -field  vector,  are  presumed  to  be  random  in  the 
experiment  described  by  Helliwell  and  Katsufrakis,  but  are  organized  to  an 
unknown  degree  in  the  phase-reversal  experiment. 

Although  there  is  not  general  agreement  that  the  modulation  period 
is  a direct  measure  of  the  phase  bunching  period  [Nunn,  1975;  Helliwell 
and  Crystal,  1 975j  we  will  assume  that  it  is  and  take  the  radian  frequency 
of  the  modulation  for  the  phase  bunching  frequency,  w = 2w/T  = 245  sec”*. 
Then,  the  value  for  a obtained  from  the  TVLF  transmission  experiment  is 


0.12. 

Again,  one  must  be  cautious  in  comparing  the  results  of  the  TVLF 
measurements  with  the  computer  simulation.  The  characteristic  plasma 


parameters,  listed  in  Table  1,  differ  considerably.  Also,  some  degree 
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Table  1 


Comparison  of  parameters  (wp,  plasma  frequency,  and  «c» 
cyclotron  frequency)  in  the  computer  simulation  study  and 
the  TVLF  phase  reversal  experiments. 


of  phase  coherence  existed  at  the  time  the  phase  reversals  were  applied 
during  the  TVLF  experiment.  Hence,  the  initial  conditions  in  the  two 
experiments  were  not  alike.  To  be  strictly  correct,  a simulation  is  required 
which  includes  periodic  phase  reversals  for  waves  propagating  in  a plasma 
in  a non-uniform  magnetic  field. 

The  magnetic  field  intensity  of  the  whistler-mode  wave  at  saturation 
may  be  estimated  if  we  identify  the  modulation  period  obtained  from  the  data 
from  the  present  experiment  with  the  phase  bunching  period  T given  by 
[He  lli  we  11,  1967] 


where  m is  the  mass  of  the  electron,  X.  is  the  wavelength  of  the  whistler- 
mode wave,  q is  the  charge  of  the  electron,  vN  is  the  velocity  of  the 
resonant  electrons  parallel  to  the  geomagnetic  field,  and  a is  the  pitch 
angle  of  the  electron  velocity  vector.  It  is  assumed  that  the  wave  is 
propagating  parallel  to  the  geomagnetic  field. 

The  plasma  density  required  to  evaluate  K and  v(|  can  be  obtained 
from  single-and  multiple-hop,  lightning -gene rated  whistlers  which  were 
occurring  throughout  this  time  period.  Propagation  through  only  one  duct 
is  apparent  in  the  data.  The  nose  frequency  obtained  by  the  nose  extension 
method  was  11.1  kHz  and  the  time  delay  at  the  nose  was  1.05  sec. 

This  places  the  duct  at  L = 3.  1 for  the  gyrof  requency  electron  den- 
sity model  [Helliwell,  1965],  The  travel  time  at  6.6  kHz  would  be  1.  14  sec 
in  good  agreement  with  the  delay  observed  for  the  arrival  of  the  phase  shifts 
at  Dunedin,  namely  1.  15  sec.  The  second  marks  along  the  axis  of  abscissas  in 
Fig.  3 are  properly  positioned  at  the  Universal  Time  seconds  in  relation  to  the  data. 
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The  minimum  amplitude  attained  during  the  phase  reversal  immediately 
following  the  mark  at  12:43:34  UT  is  measured  to  occur  at  12:43:34.  15  UT. 
This  phase  reversal  was  transmitted  ai  12:43:33  UT.  The  electron  density 

at  the  equator  obtained  from  the  travel  time  of  the  whistlers  at  the  nose 

-3 

requency  is  345  cm  . The  index  of  refraction  for  a 6.6  kHz  wave  prop- 
agating parallel  to  the  geomagnetic  field  is  then  13*8  which  yields  a wave- 
length of  3290  m.  For  these  plasma  parameters  the  parallel  velocity,  v,|  , 

ox  an  electron  which  is  resonant  with  the  wave  through  a cyclotron  resonance 
7 

is  7.  3 X 10  m/sec.  Helliwell  (_  1 96 7 J has  shown  that  the  characteristic 
pitch  angle  for  the  interaction  is  30  . Since  the  measurements  do  not  provide 
this  pitch  angle,  30  will  be  adopted  for  these  calculations.  Substituting 
f°r  T,  v and  a in  (6),  the  value  obtained  for  in  the  interaction  region 
is  3.  5 pT  (mV). 
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CONCLUSIONS 

We  have  presented  the  results  of  a controlled  VLF  transmission 
experiment  conducted  to  determine  the  role  played  by  the  phase  of  a 
whistler-mode  wave  in  wave-particle  interactions. 

It  is  observed  that  the  phase  of  the  signal  remains  coherent  through- 
out the  magnetospheric  propagation  path  to  the  receiver  in  the  conjugate 
region  of  the  transmitter. 

When  periodic  phase  reversals  are  applied  to  the  transmitted  signal, 
the  response  of  the  magnetosphere  is  analogous  to  that  of  a tuned  electric 
circuit  with  a bandwidth  comparable  to  bandwidths  predicted  for  cyclotron 
resonance. 


In  this  experiment  the  amplitude  saturated  at  a magnetic  field 


intensity  of  3.5  pT  (mY),  Following  a phase  reversal,  the  "linear"  growth 


rate  is  found  to  be  0.12  times  the  phase  bunching  frequency  at  saturation 
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THE  IVAN  A.  GETTING  LABORATORIES 


The  Laboratory  Operation*  of  The  Aeroapace  Corporation  i*  conducting 
experimental  and  theoretical  investigation*  necessary  for  the  evaluation  and 
application  of  ecientific  advances  to  new  military  concept*  and  system* . Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos - 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronic*  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Material*  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanic*  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Science*  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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